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Abstract The limit theorems of discrete- and continuous-time quantum walks 
qj , on the line have been intensively studied. We show a relation among limit 

distributions of quantum walks, Heun differential equations and Gauss dif- 
ferential equations. Indeed, we derive the second-order Fucksian differential 
equations which limit density functions of quantum walks satisfy. Moreover, 
using both differential equations, we discuss a relationship between discrete- 
and continuous-time quantum walks. Taking suitable limit, we can transform 
a Heun equation obtained from the limit density function of the discrete-time 
quantum walk to a Gauss equation given by that of the continuous-time quan- 
tum walk. 
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1 Introduction 

The discrete-time quantum walk (QW), which is a quantum counterpart of 
the classical random walk, has been extensively investigated since Ambainis 
et al. [5] studied a detail of the walk. The continuous-time QW was proposed 
■ by Farhi and Gutmann [5], and has been analyzed on not only regular graphs 

but also complex networks (e.g. [SHH1I32)- A further development in the theory 
of the QW during recent 10 years showed us several novel properties on both 
QWs. The behavior of the QW is quite different from that of the random 
walk and is expected to be connected with various phenomena in quantum 
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mechanics. A relation between the QW and quantum computer has been also 
discussed (e.g. [IJ[2) ( 23 ). The quantum search algorithm designed by the QW 
is one of the important applications. The Grover search algorithm can be 
considered as a discrete-time QW on the complete graph and it produces 
speed-up at the square-root rate for the corresponding classical search. There 
are some reviews of the QW [TUllTTlfrSllgT] ■ 

One of the purposes in the study of the QW is to derive the limit dis- 
tribution and the asymptotic behavior as time tends to infinity. The limit 
distribution has been obtained for various kinds of QWs. In the present paper, 
we concentrate on the QW in the one-dimensional space. Observing previous 
results, one can mention that each limit distribution usually has a compact 
support and admits a singularity at the boundary of the support. This is an 
interesting property of the limit distribution, because the QW does not pos- 
sess any singularity in space. To the fact, we now address a question: How to 
understand the singularity of the QW? Motivated by this question, we inves- 
tigate the limit distribution of the QW from a viewpoint of the differential 
equation. 

We begin with finding a differential equation for the limit distribution. 
For example, consider discrete- and continuous-time symmetric simple random 
walks on the line. The central limit theorem for the walk, that is, lim^oo P(Y t /^/t < 
x) = f_ e~ v I 2 /V2tt dy, is well-known, where Y t denotes the walker's posi- 
tion at time t and V(Y t — x) is the probability that the walker is at position 
x at time t. For the limit density function fix) = e~~ x ^ 2 /\/2n, an equation 
d f(x)/dx + xdf(x)/dx 4- f(x) — is derived. This equation comes from 
the theory of the diffusion equation du(x,t)/dt = d 2 u(x,t)/dx 2 , and in par- 
ticular, it is also obtained by the self- similarity of the fundamental solution 
u(x,t) = [1/Vi)f(x/Vi). 

Let us consider both discrete-time and continuous-time QWs on the line. 
The classical random walk corresponds to the diffusion process, while for the 
QW the corresponding process is not known. Thus, to find out the process 
would remain as one of important problems. Taking an account of the singu- 
larity in the limit distribution, we treat a class of Fucksian linear differential 
equations of the second order. For discrete-time (resp. continuous-time) QWs, 
we are led to a Heun's differential equation (HE) (resp. a hypergeometric dif- 
ferential equation by Gauss (the Gauss equation, GE)). The GE is one of 
typical Fucksian equations and represents the Fucksian equations with exactly 
three regular singular points. On the other hand, the HE was proposed in 1888 
by Heun [7] and admits four regular singular points. 

Moreover, a concept of confluence in between the HE and the GE helps 
us to understand a relationship between discrete- and continuous-time QWs 
clearly. Through these results, authors believe that the HE and the GE play 
an important role in understanding the QWs. 

To this end, we will mention a significant remark on the QW and the 
BC\ Inozemtsev model, which is an integrable quantum system. These two 
models are to be connected through HEs. Therefore, we would like to expect 
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a new mathematical theory on the relation between QWs, HEs, and quantum 
mechanics. 

The present paper is organized as follows. In Sect. [2l we explain both the 
one-dimensional discrete-time QW and the HE. After introducing the limit 
distribution obtained by Konno [T2UT1] . we discuss a relation between the 
discrete-time QW and the HE. In Sect. [31 we concentrate on the continuous- 
time QW, and show that the GE relates with the limit density function of 
the walk. Section [4] is devoted to a connection between the HE and the GE 
in order to understand a relationship between discrete- and continuous-time 
QWs given by Strauch [24] , In the final section, we summarize the relations 
obtained here and propose a future problem on our results. Furthermore, we 
transform the HE to the BC\ Inozemtsev system in |Appcndix 1 and show 
relations among measure of the QW, the HE and the GE in Appendix 2| 



2 A relation between the discrete-time QW and the HE 



In this section, we discuss a relation between the discrete-time QW in one 
dimension and the HE. At first we define the QW on the line. Let \x) (ieZ = 
{0, ±1, ±2, . . .}) be infinite components vectors which denote the position of 
the walker. Here, x-th component of |x) is 1 and the other is 0. Let \i/jt(x)) € C 2 
be the amplitude of the walker at position x at time t £ {0, 1,2,.. .}, where C 
is the set of complex numbers. The amplitude of the walk at time t is expressed 

by _ 

|# t > = 5>>® (i) 

The time evolution of the walk can be defined by the following unitary matrix: 



U = 



a b 
c d 



where a, b,c,d € C. Moreover, we introduce two matrices: 

P 

Note that P + Q = U . Then the evolution is determined by 
\ijj t +i{x)) = P |Vt(z + 1)> + Q \M* ~ 1)> • 



a b 




'0 0" 





,Q = 


c d 



(2) 



(3) 



(4) 



The probability that the discrete-time quantum walker x[ d ^ is at position x 
at time t, ¥(xj: d ' > = x), is defined by 



V(X ( t d) 



x) = (ip t (x)\ipt(x)) , 



(5) 



where (tpt(x)\ denotes the conjugate transposed vector of \ipt{x)). In the present 
paper, we take the initial state as 



\Mx)) 



T [l/V2, i/y/2] (x = 0), 
T [0,0] (x^0), 



(G) 
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where T is the transposed operator. Equation © is well-known as the initial 
state that gives the symmetric probability distribution about the origin [15 . 

For the walk, some limit theorems and asymptotic behaviors have been 
obtained. In particular, we focus on the density function of probability distri- 
bution as t — > oo. The limit distribution of the QW was given by Konno (T2J 
03]. That is, for abed ^ 0, we have 

/ x (d) \ fX 
lim P M- < x )= / f {d) (y)l { -\au\a\)(y)dy, (7) 



t— >OQ 



where 



/«>(*)= z 4=S =^. (8) 



7r(l-^)v1 
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and Ia(x) = 1 if x G A, /^(x) = if a; <^ A. A lot of limit distribu- 
tions of discrete-time QWs are often described by using this function f^ d '(x) 
(e.g. pfl CTTMTTl[TBll!Z2] ) . Figure [1] depicts the comparison between the proba- 
bility distribution at time 500 and the limit density function for the walk with 
a = b = c = —d = l/y/2, which is called the Hadamard walk. 




Fig. 1 Comparison between the probability distribution and the limit density function of 
the discrete-time QW with a = b = c = —d = l/\/2. Figure (a) is the probability distribution 
at time t = 500. Figure (b) is the limit density function of the walk. 



From now, we discuss a relation between the limit density function f( d \x) 
and the HE. The canonical form of the HE is given by 

d 2 u f 7 S e \ du a/3z — q 

d^ + \~z + V^l + J~6) d^ + z(z- l)(z °' (9) 

where 9 G C is one of the singular points and a, ft, 7, 5, e, q € C. The five 
parameters a, /3, 7, S, e are linked by the relation a + /3 + l = 7 + 6 + e and 
the parameter q is called the accessory parameter. The HE has four regular 
singularities at z — 0, 1,9, 00. Particularly the HE with j = S = e = l/2is 
called Lame's equation and it have been analyzed in detail. In some special 
cases, the HE becomes the GE (e.g. 9 = l,q= a/3 case). 
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We find that f^{x) satisfies the following differential equation: 

(1 _ ^)(| a | 2 _ x 2 )^ d \x) - x(4|a| 2 + 3 - lx 2 )j-f d \x) 

+(9a; 2 -2|a| 2 -l)/W(a;)=0. (10) 

By a change of an independent variable x 2 — z, we obtain one of our main 

results: 

Theorem 1 

^„« (z)+ (i + _i_ + _i_)|„« (z) 

£,_ 2|a| 2 + l 

uW(«)=0, (11) 



4 Z 4 



z{z-l){z-\a\>) 



whe 



M W = ^ /1 -' a ' 2 . (12) 

7r(l-2!)v1HF^ 

We should remark that Eq. (jlll) is equivalent to the HE with 

3 1 3 2|a| 2 + l l9 , . 

a = P = -, 7 = S = 2, e = -, q = -LL , Q = \ a \ 2 . ( i3 ) 



3 A relation between the continuous-time QW and the GE 



z(z-l) — + {{a + 13 + \)z- 7}—+ a/3u = (14) 



In this section we will show a relationship between the continuous-time QW 
on Z and the GE: 

, d 2 u n „. , du 

,_ + { (a + /> + l),- 7 }- 

where a, /3, 7 £ C are parameters. Equation (|14[) is a second-order Fucksian 
equation with three regular singularities at z = 0, 1, 00. 

We give a definition of the continuous-time QW on the line. At first, we 
consider the amplitude ipt{x) G C at position x £ Z at time t(> 0) instead 
of ^(a;)) £ C 2 for the discrete-time QW. The evolution of the amplitude is 
defined by 

where v > 0. The probability that the walker is at position x at time t is 
denoted by P(X^ = x) = \ipt{x)\ 2 , where X\°^ is the continuous-time quan- 
tum walker's position at time t. We take an initial state as ipo (0) = 1 and 
ipo(x) — (x 7^ 0). Konno [T3] and Gottlieb [5] got the limit theorem for the 
continuous-time QW as follows: 



hm P M- < x = / / (c) (y) J ( _ 2 ,, 2 „) (») (16) 
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where 

f (c) (*)= , * 2 - (17) 

In Fig. [21 we show the comparison between the probability distribution at time 
500 and the limit density function for the walk with v = l/2y/2. 




Fig. 2 Comparison between the probability distribution and the limit density function of 
the continuous-time QW with u = l/2\/2. Figure (a) is the probability distribution at time 
t = 500. Figure (b) is the limit density function of the walk. 



We see that the function f^ c \x) satisfies the following differential equation: 

(4^ 2 - * 2 )^/ (c) (x) - 3z^/ (c) (z) ~ &\x) = 0. (18) 

By the change of an independent variable x 2 /4u 2 = z, we have a continuous- 
time counterpart of Theorem [TJ 

Theorem 2 

z(z - l)^ C \z) + (<2z ±) j-/ c Hz) + \u^{z) = 0, (19) 

where 

u^\z) = - )==. (20) 

Note that Eq. is the GE with a = /3 = 7 = 1/2. 



4 A relation between discrete-time and continuous-time QWs in 
the theory of the differential equation 

In this section, we discuss a relation between discrete- and continuous-time 
QWs on the second-order differential equations which were obtained in previ- 
ous sections. Strauch |24] found a relation between both walks by transforming 
Eq. flU to Eq. (fT5]l in a suitable limit. D'Alessandro [I] focused on both QWs 
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on general graphs and derived the dynamics of continuous-time walks as a 
limit of discrete-time dynamics. 

Considering a confluent type for the HE which the function vS d > (z) satisfies, 
we can also obtain the GE which the function u^ c \z) satisfies. For a scaling 
parameter r, substituting z = t/r in Eq. (jlip . we have 

At{\a\\ -t)(l- 1) ^v(t) + 2 { ^ - (5|a| 2 + 4)t + |a| 2 r| ±v(t) 

+ (~ - 2|a| a - = 0, (21) 

with u(t) = vS d '{t/r). As t -> oo, |a| — > under the condition \a\ 2 r — > 1, we 
obtain a confluent HE: 

t(t-D-^ + (2i--j- + r = 0. (22) 

Equation (j2"2")l is equivalent to the GE which was obtained from the continuous- 
time QW (see Eq. ((T5jl ). This result corresponds to the result in Strauch |24) . 
We can confirm his result in a relation between the HE and the GE. 



5 Summary 

In this section, we discuss and conclude our results. In the present paper, 
we found that for the discrete-time QW, the limit density function f^ d \x) 
satisfies a HE, while for the continuous-time QW, the limit density function 
f( c \x) satisfies a GE. Moreover, considering the confluent HE, we confirmed 
a relation between the discrete- and the continuous-time QWs corresponding 
to the result by Strauch |2~4"] . 

As a significant remark, we would like to mention a relationship between the 
discrete-time QW and a Schrodinger equation through the HE (see |Appendix 1| 
for detail). The Hamiltonian of the continuous-time QW is given by the adja- 
cency matrix of the graph on which the walk is defined [5] . However the Hamil- 
tonian of the discrete-time QW is not known. On the other hand, it was shown 
by Takemura {2"51l26U27U2gll29U30| that the HE can be transformed to the BC X 
Inozemtsev system which is a one-particle integrable quantum system. In ad- 
dition, the BCi Inozemtsev system includes the Calogero-Moser-Sutherland 
system or the Olshanetsky-Perelonov system [21] . And the Hamiltonian is ex- 
pressed with the Weierstrass p-function. Thus, to find a relation among the 
discrete-time QW, the HE and the BC\ Inozemtsev system might be one of 
the interesting problems. 



Appendix 1 A relation between the discrete-time QW and the 
BC± Inozemtsev model 

In this Appendix, wc transform the HE obtained in Sect. [2] to the BC\ Inozemtsev system. 
The BCi Inozemtsev model is known as an integrable one-particle quantum system including 
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the Calogero-Moser-Sutherland system (see 25 26 27,28 29,30.). The Hamiltonian H of the 
BC± Inozemtsev model is defined by 

H = --^+Y J {l J +l)p(x + w J ), (23) 

i — 

where the function p{x) is the Weierstrass p-function with periods (2wi , 2W3) and the 
parameters lj 6C(j = 0, 1, 2, 3) are constants. We should note that Wi, 103 £ C are linearly 
independent on R, where M is the set of real numbers. By using the following transformation 
for Eq. QTJ (see Takemura [27ll28l for detail), 



(24) 



z _ p(x) - p(wi) 
p(ma) - 

and putting g(x) = — lal 2 ) 3 / 4 , we can get 

ff<?0) = 2 ~^ 2 {pN - p(-wi - w 3 )} 9(x), (25) 

where 

d 2 1 3 

-ff = _ T^O") + jP( x — u>i — 103) + 2p(x + w 3 ). (26) 

The operator H is the Hamiltonian of the BC\ Inozemtsev model with ?o = — ■k > h = 0, (2 = 
— f 1 ^3 = 1) ^0 = 0, u>2 = — lUi — U13, where uii , U13 are arbitrary non-zero complex numbers 
and linearly independent on R. Therefore the eigenvalue of the Hamiltonian H determined 
by Eq. i !2!)[ i is — ^LJ_ |p(to 1 ) — p(—wi — U13)} and the eigenfunction corresponding to the 
eigenvalue is the above mentioned function g(x). Via the HE, we found that the limit density 
function (x) of the discrete-time QWs is related to the BC\ Inozemtsev model. 



Appendix 2 Relations among measure of the QW, the HE and the 
GE 

In this Appendix, we focus on both limit measures / ( d ) (x) dx and /( c ) (x) dx and discuss a 
relation among the QW, the HE and the GE. Putting x 2 = z for the measure /< d ) (a;) dx of 
the discrete-time QW, we get 



(X)dx : 



2n(l-z)y / \a\ 2 z- 



27r(l-«)VW 



■-dz (x > 0), 
= dz (x < 0). 



(27) 



The function w^iz) = y'l - \a\ 2 /2ir(l - z)^/\a\ 2 z - z 2 satisfies the HE with a = (5 



2:7=1 
d 2 



8 = 2, 



3|n|-+2 



M 5 



dz 2 



w {i) (z) + I 2 +^— + 



z z — 1 z — \a\ 2 I dz 



. 3|o| 2 +2 
— w {d \z) -\ 2 



z{z-\){z-\a\*) 



w < - d \z)=0. (28) 



In the case of the continuous-time QW, putting x 2 /4u 2 = z for f^(x) dx, we obtain 



/< c > (x)dx 



0= dz (x > 0), 

- _ z 2 

dz (x < 0). 



7T\/ 2 — 2 



(29) 



7TV 2-2 4 
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For the function w^ c \z) = v/nV z — z 2 , the following GE with a = /3 = 1,7 = % is realized: 

Z(z- l)^-w^ c \z) + (zz- -~) —w^{z)+w ( - c \z) = 0. (30) 
dz A \ 2 J dz 

Note that we can derive Eq. (1306 from Eq. (1 2S 6 in a similar fashion as in Sect. [4] More- 
over Eq. (1281) can be transformed into the eigen equation of the BC\ Inozemtsev model 
with Iq = —^,h = —1,1.2 = ~§i^3 = ~^, w = 0,t02 = — tt>i — u>3 and the eigenvalue 

2— lal 2 

12 IPW - p(-u>i - W 3 )}. 
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